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Abstract 

We compute the masses and decay widths of the gluonia using QCD spectral sum rules and low-energy theorems. 
In the scalar sector, one finds a gluonium having a mass Mc = (1-5 ± 0.2) GeV, which decays mainly into the 
U{1)a channels r/77' and 47r'^. However, for a consistency of the whole approach, one needs broad-low mass 
gluonia (the cjb and its radial excitation), which couple strongly to the quark degrees of freedom similarly to 
the rj' of the U{1)a sector. Combining these results with the ones for the qq quarkonia, we present maximal 
gluonium-quarkonium mixing schemes, which can provide quite a good description of the complex spectra and 
various decay widths of the observed scalar mesons (t(1.), /o(0.98), /o(1.37), /o(1. 5) and /,/(!. 71). In the tensor 
sector, the gluonium mass is found to be Mt — (2.0 ± 0.1) GeV, which makes the C(2.2) a good 2++ gluonium 
candidate, even though we expect a rich population of 2+"*" gluonia in this region. In the pseudoscalar channel, 
the gluonium mass is found to be Mp ~ (2.05 ±0.19) GeV, while we also show that the E/l(1A4) couples more 
weakly to the gluonic current than the 77' (0.96), which can favour its interpretation as the first radial excitation 
of the 77'(0.96). 
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1 Introduction 



There is now clear evidence from many processes that QCD is the theory of strong interactions. All hard 
processes satisfy the asymptotic freedom property of QCD, while the complicated hadron properties can be 
explained by different non-perturbative methods (effective Lagrangian, lattice calculations, QCD spectral sum 
rule (QSSR),...). However, in addition to the well-known mesons and baryons bound states, one of the main 
consequences of the QCD theory is the possible existence of the gluon bound states (gluonia or glueballs) 
or/and of a gluon continuum. The theoretical interest in the gluonia sector has started a long time ago, since 
the pionneering work of Fritzsch and Cell-Mann [Q, as shown by the long list of publications on this topic. 
However, despite these different efforts, the theoretical and experimental status remains unclear. From the 
theoretical point of view, this is due to our poor control of the gluon dynamics. In particular, there is not 
yet any convincing dynamical approach that computes the mixing of the gluonia with quarkonium states (see 
however [|[ jsj), which is necessary if one wishes to make contact with the experimental data. At this level, 
only phenomenological scheme is available, where the mixing angle is only fitted from the data. From the 
experimental point of view, the difficulty also arises in the same way, as the observed resonances can be a 
mixing between gluonia and quarkonia. However, some selected processes such as the J /tp radiative decays can 
favour more the production of gluonia than of quarkonia, while the measurement of the two-photon widths has 
been used for a long time as a good gluonia signature. However, as proposed some years ago in Q, a good 
signature for the presence of the gluonia can be obtained from the ratio of the previous two quantities referred 
to as "stickiness" . Another possibility for signing the nature of an almost pure gluonium is the measurements 
of its width into the U(l)-like channels: 77'??', jy??', and 477", which are expected to be more dominant than, for 
instance, the one into pair of pions, leading [|[ ^ to the conclusions that the C(1.6) observed by the GAMS 
group is an almost pure glueball state. The different experimental progress done during these last few years 
[^, though not yet very conclusive , is encouraging for improving the theoretical predictions. Some recent 
efforts have been accomplished in this direction from lattice calculations However, the "apparent" 

disagreements of these lattice results can simply reflect the true systematic error of the estimates from this 
approach. Recently, some QCD inequalities among the gluonia masses have been derived 12 1. Some qualitative 
phenomenological attempts | |l3t have also been proposed for explaining the nature of the scalar meson seen by 
the Crystal Barrel collaboration Motivated by these different steps forward, we plan to give in this paper 
an almost complete scheme and an update of the predictions from QSSR a la SVZ [ p^ ^, where, in my opinion, 
the results have often been misinterpreted in the literature, and in some cases ignored. A long list of sum rules 
results exists in the literature [0-||2^, which needs to be updated because of the progress accomplished in 
QCD during the last few years. This is the aim of this work. The paper is organized as follows: 
In section 2, we give a general discussion on the gluonia and the classification of different currents. In section 
3, we give a short introduction to the method of QSSR. In sections 4 and 5, we update mainly the work of |^ 
(hereafter referred to as NV) for the 0++ sector by improving the QSSR predictions on the masses and couplings 
of the gluonia, and by extending the uses of some low-energy theorems (LET) for their decay widths, into new 
predictions for their decay into 47r°. In section 6, we compute the properties (masses and widths) of the scalar 
quarkonium. In section 7, as in p7| (hereafter referred to as BN), we use some maximal quarkonium-gluonium 
mixing schemes for explaining the scalar mesons data below and above 1 CeV. In section 8, we discuss the QSSR 
estimate of the 2++ gluonium mass and coupling, which is an update of the work of ||l8| (hereafter referred to 
as SN). We close this section by giving the widths of the 2++ gluonium using a meson-gluonium mixing scheme 
and LET as in |]. 

In section 9, we discuss the QSSR estimate of the gluonium mass and coupling (update of (l8|). We use 
the quarkonium-gluonium mixing scheme for predicting the 77 and p7 decay of the pseudoscalar gluonium. An 
attempt to explain the property of the E / 1, (1.44) is given. A summary of our results is given in the two tables 
in section 10. 



2 The gluonic currents 

In this paper, we shall consider the lowest-dimension gluonic currents that can be built from the gluon fields 
and are gauge- invariant: 



Js = /3(a,)G„0G"^, 

^For a recent review on the sum rules, see e.g. |l6t 
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= (^)trG„^G'"^ (1) 
where the sum over colour is understood and in our notations |l6| ^: 



Gpjy — -:^^t_LuapG°'^ . (2) 

These currents have respectively the quantum numbers of the J^^ = 0++, 2++ and gluonia ^, which 
are famihar in QCD. The former two enter into the QCD energy- momentum tensor 6^^, while the later is the 
11(1)^1 axial-anomaly current. The renormalizations and corresponding anomalous dimension of the scalar and 
pseudoscalar currents have been studied in QCD where renormalization group invariant quantities have 
been built. In the approximations, at which we are working, we can ignore renormalization effects relevant at 
higher orders of the as radiative corrections. Besides the sum rules analysis of the corresponding two-point 
correlators, we shall also study the gluonia properties using some LET, while more phenomenological mixing 
schemes will be presented to explain the data. 



3 QCD spectral sum rules 

The analysis of the gluonia masses and couplings will be done using the method of QSSR. In so doing, we shall 
work with the generic two-point correlator: 

V'G(g') = ijd^x e'«^ (0|rjG(x) (Jg(0))^ |0), (3) 

built from the previous gluonic currents Jg{x). Thanks to its analyticity property, the correlator obeys the 
well-known Kallen-Lehmann dispersion relation: 

^G{q^)= ^ - Imi^ait) + (4) 

where ... represent subtraction points, which are polynomials in the g^-variable. This sum rule expresses in 
a clear way the duality between the integral involving the spectral function lu\il>Q[t) (which can be measured 
experimentally), and the full correlator il)Q{q^), which can be calculated directly in QCD using perturbation 
theory and the Wilson expansion, provided that —q^ is much greater than A^. 



3.1 The Operator Product Expansion 

By adding to the usual perturbative expression of the correlator, the non-perturbative contributions, as parametrized 
by the vacuum condensates of higher and higher dimensions in the OPE ||l5|| , the two-point correlator reads in 
QCD: 

^G{q^)^ E , L/2 - E C(^)(g^z.)(OM), (5) 

D=0,2,4,... ( 1 ) dimO=D 

where v is an arbitrary scale that separates the long- and short-distance dynamics; C^"^^ are the Wilson coef- 
ficients calculable in perturbative QCD by means of Feynman diagrams techniques. In the present analysis, 
we shall limit ourselves to the computation of the gluonia masses in the massless quark limit ^, which we shall 
compare with the results obtained from other methods such as lattice calculations and QCD inequalities. The 

^The corresponding differential equation for the running coupling is dds/dt = as/3(«s) where t = 1/2 log(— g^/z/^). 

^We shall not consider the pseudotensor 2""*" in this paper. A QSSR analysis of this channel is in progress. 
Quark mass corrections which will appear at higher order through internal fermion loops in the OPE at the order 
a g {fh g / ) log{— q-^ /u^), are expected to be tiny at the scale (see next sections) where the sum rules are optimized {fhs being 
the running quark mass). 
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OPE is shown diagramatically in Fig. 1: 





+ . 



Figure 1: Diagrammatic representation of the OPE of the gluonic two-point correlator 



• The D — operator corresponds to the usual case of the nai've perturbative contribution. In the massless 
light quark limit , some additionnal D = 2 terms not included in the OPE may eventually manifest after 
the resummation of the QCD perturbative series However, due to the unclear present status of 
this contribution, and to the inaccurate quantitative control of the contribution of this term from ultraviolet 
renormalons calculus, we shall not explicitly consider this effect in our analysis. Alternatively, we estimate, 
like in the case of r decay ^3| - ||3^ , the next unkown higher-order correction using a geometrical growth of 
the QCD series, which value is in agreement with the one from experimental measurements |36 and from 
an estimate |]37| based on a i?S'-invariant approach ^ within the principle of minimal sensitivity |38|. Then, 
we argue, like in the case of r decay, that the errors induced by the truncation of the QCD series are given 
by this last term (theorem of divergent series [^), which we have multiplied by a factor 2 as in |35| for a 
conservative estimate. The errors estimated in this way are comparable in strength with some of the previous 
UV renormalon estimates, and should already include in it this rcnormalon contribution (manifestation of the 
resummed QCD series). However, as mentioned previously, the UV renormalon estimate should be considered 
to be very qualitative at the present status of the art (dominance of the second chain of renormalons js^,...). 

• The D = A dimension operators, which will give the dominant contributions in the chiral limit m.^ 
is the gluonic condensate (ofsG^), introduced by SVZ and which has been estimated recently from the 
e~^e~ — > / = 1 hadron data | |4C| | and from the heavy quarks mass sphttings [HI: 



(0.07 ±0.01) GeV^ 



(6) 



in agreement with different post-SVZ estimates quoted in these papers. 

• The D = 6 contribution is dominated by the triple gluon condensate gfabc{G°'G^G'^), whose direct extraction 
from the data is still lacking. However, one can estimate its approximate value from the dilute gas instanton 
model 

g^'fabciG'^G'G^) « (1.5 ± 0.5)GeV2)(a,G2) , (7) 
where we have given the error in such a way that the estimate is valid within a factor 2. 

• The validity of the SVZ expansion has been intensively studied in the theoretical literature while the 
good agreement for the values of the QCD coupling as from r decay |33 and LEP data f\, the observation 
of the running of as at low-scale p3|, and the agreement of the measured QCD condensates from semi-inclusive 
tau decays and spectral moments |44 with the ones from phenomenological fits can be considered as 
its phenomenological confirmation. 

• In the case of gluonia studied in this paper, it has also been argued jl^, using the dilute gas approximation, 
that in some channels, instanton plus anti-instanton effects manifest themselves as higher dimension (D = 11) 
operators. On the other hand, its quantitative estimate is quite inaccurate because of the great sensitivity of 

®The iiS-invariant approach can be applied in the present case of gluonia as the corresponding two-point function starts to 
lowest order as , such that the unknown correction to be estimated is at order . 

'^New data from deep inelastic scattering give now a value of a,s in a better agreement with the two former |42|. 
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the result on the QCD scale A, and on some other less controllable parameters and coefBcients. However, at 
the scale (gluonia scale) where the following sum rules are optimized, which is much higher than the usual case 
of the p meson, we consider that one can safely omit these terms like any other highcr-dimcnsional operators 
beyond Z) = 8. 

3.2 The spectral function and its experimental measurement 

The experimental measurement of the spectral function is best illustrated in the case of the flavour-diagonal 
light quark vector current, where the spectral function Imn(t) can be related to the e'^e^ into 7=1 hadrons 
data via the optical theorem as: 

cr(e+e~ ^ / = 1 hadrons) = - Imn(t). (8) 

t TT 

One can also relate the spectral function to the leptonic width of the p resonance: 

via the meson coupling to the electromagnetic current: 

W\P) = (10) 

within a vector meson dominance assumption. More generally, one can introduce the decay constant /g analo- 
gous to = 93.3 MeV: 

(0|Jg|G) = V2/gA4 (11) 
where ... represent the Lorentz structure of the matrix elements. 

3.3 The form of the sum rules 

On the phenomenological side, the improvement of the previous dispersion relation comes from the uses of an 
infinite number of derivatives and infinite values of g^, but keeping their ratio fixed as r = n/g^. In this way, 
one obtains the Laplace or Borel or exponential sum rules [|l^ ^ : 

Cg{,t) = / dt exp(-ir) - ImV'G(i), (12) 

where t< is the hadronic threshold. The advantage of this sum rule with respect to the previous dispersion 
relation is the presence of the exponential weight factor, which enhances the contribution of the lowest resonance 
and low-energy region accessible experimentally. For the QCD side, this procedure has eliminated the ambiguity 
carried by subtraction constants, arbitrary polynomial in g^, and has improved the convergence of the OPE by 
the presence of the factorial dumping factor for each condensates of given dimensions. 
The ratio of sum rules: 

7^G = -|-log/:G, (13) 

or its slight modification, is a useful quantity to work with, in the determination of the resonance mass, as it is 
equal to the mass squared, in the simple duality ansatz parametrization: 

"one resonance" 5 - M^) + "QCD continuum" 6 (i - tc), (14) 

of the spectral function, where the resonance enters by its coupling to the quark current; tc is the continuum 
threshold which is, like the sum rule variable r, an a priori arbitrary parameter. As tested in the meson channels, 
this parametrization gives a good description of the spectral integral, in the sum rule analysis. In some cases, 
we shall also use finite energy sum Rule (FESR) Q 0: 



A^[?) ^ r dtt-- lu^i^ait), ^ (15) 



((») 

where n is an integer, in order to double check the estimate obtained from the Laplace sum rules. 
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3.4 Conservative optimization criteria 

Different optimization criteria are proposed in the literature, which, to my opinion, complete one another, if 
used carefully. The sum rule window of SVZ is a compromise region where, at the same time, the OPE makes 
sense while the spectral integral is still dominated by the lowest resonance. This is indeed satisfied when the 
Laplace sum rule presents a minimum in r, where there is an equilibrium between the non-perturbative and 
high-energy region effects. However, this criterion is not yet sufficient as the value of this minimum in r can 
still be greatly affected by the value of the continuum threshold <c- The needed extra condition is to find the 
region where the result has also a minimal sensitivity on the change of the tc values {tc stability) . The tc values 
obtained in this way are about the same as the one from the so-called heat evolution test of the local duality 



FESR 1 47 . However, in some cases, this tc value is too high, compared with the mass of the observed radial 
excitation, and the procedure tends to overestimate the predictions. More precisely, the result obtained in this 
way can be considered as a phenomenological upper limit. Therefore, in order to have a conservative prediction 
from the sum rules method, one can consider the value of tc at which one starts to have a r-stability up to 
where one has a tc stability. In case there is no tc stability nor FESR constraint on tc, one can consider that 
the prediction is still unreliable. In this paper, we shall limit ourselves to extracting the results satisfying the 
T (Laplace) or n (FESR) and tc stability criteria.^ 

4 Mass and decay constant of the 0++ scalar gluonia 

4.1 The gluonium two-point correlator in QCD 

We shall be concerned with the correlator: 

^M^) ^ m I d^x e^-^- {0\Te>^Jx)em^\Q), (16) 



where d^n is the improved QCD energy-momentum tensor (neglecting heavy quarks) whose anomalous trace 
reads, in standard notations: 

e'^ix) = i/3(a,)G2-f (l+7„(a3)) ^ m,V;,^,. (17) 

u,d.s 

Its leading-order perturbative and non-perturbative expressions in as have been obtained by the authors of 
[ p7[ . To two- loop accuracy in the MS scheme, its perturbative expression has been obtained by while the 
radiative correction to the gluon condensate has been derived in | |2^ . Using a simplified version of the notation 
in Eq. (5): 

lPs{q^)^ J2 Cd{Od), (18) 
_D=0,4,... 

one obtains for three flavours and by normalizing the result with (/3(as)/as)^: 

2 



Ci{Oi) 



C6(06> - 2(^)|l-^a.log-^|5Vabc(G''G''G^) 

G8(08) = U{{asfabcGl^Glpf)-{{aJabcG%Gl^f). (19) 
From its asymptotic behaviour {q^ log), we can write a twice-subtracted dispersion relation for ips{q'^): 

V'.(g') = V'.(0) + gV;(0) + 9M -—^±^—iuv^,{t), (20) 

Jo V' ~ q ^ ie) IT 

where the subtraction constant V's(O) is known from the LET to be [ p^f : 

^,(0) = -16^(a,G2). (21) 

TT 



'Many results in the literature on QCD spectral sum rules literature are obtained using only the first condition. 
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4.2 The gluonium sum rules 



Using standard QSSR technology, one can derive, from the previous expression of the two-point correlator, the 
unsubtracted Laplace sum rule (USR): 



f°° 1 

/ dt exp(— tr) — limps{t) 



^ + i ( v) (^^ + ^^^^^ " ^^^^ " " ^^s^^"^'' 



-TT— (ofsG )r +- — 5 ]abc\(j G G )t H (asG ) r 



(22) 



and the subtracted sum rule (SSR), which depends crucially on the subtraction constant V's(O): 
Ci-^^ = / - exp(-tT) - Imi^.it) 



1 + ^ (^) + 4/3i(l - je) - 4^^ log - log (tA2; 



+V'«(0)-4 



Ztt 



where: 



/?! log^^A' 

Throught this paper, we shall use for three active flavours p^ : 

A = (375 ± 125) MeV. 



(23) 



(24) 



(25) 



• One can notice that the perturbative correction is large for the Laplace sum rules £!"•' , but it tends to cancel 
in the ratio of moments. As discussed earlier, we estimate the higher-order unknown perturbative corrections 
using a geometrical growth of the QCD series. Then, we expect that the "effective" correction including 
log log terms normalized to the lowest-order perturbative graph is here about: 



Si^'^ « ±2 X 18 



V TT 



(26) 



which we consider as the error due to the truncation of the perturbative QCD series (last known term of the 
series), which should include in it some (eventual) effects of the l/Q^-term induced by the resuumation of this 
series. The factor 2 is a conservative estimate of the errors. 

• We shall see later on, that the non-perturbative contributions including the corresponding corrections are 
small at the scale where the sum rules are optimized, such that the dominant errors from the non-perturbative 
effects come from the values of the condensates. 

• One should also notice that it can be inconsistent to work here in a pure Yang-Mills, as our QCD parameters 
{as and gluon condensates) have been extracted from the data in the presence of quarks. Therefore, one should 
be careful in comparing the results in the present paper with e.g. the lattice ones. 

• The USR and the SSR ci~^'^ have been used in NV. The authors of Ref. have also used the SSR in 
order to set the sum rule scale in the scalar gluonium channel, leading to the conclusion that the scalar gluonium 
energy scale should be much larger (or r much smaller) than the one of ordinary hadrons. We shall see later 
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on that the SSR Ci satisfies this criteria, due mainly to the relative importance of the ips{Q) contribution. 
However, this conclusion is not universal, as it does not necessarily apply to the USR, where r can be slightly 
larger (or the sum rule energy scale slightly lower) than in the SSR. A simple inspection of these different 
sum rules indicates that cf'' is the only sum rule that can be sensitive to a low-mass resonance below 1 GeV, 
as it is a low moment in t and it can stabilize at large value of r. Higher moments in Cs ' ,are, on the 
contrary, sensitive to a higher mass resonance. As discussed before, ci~^^ is quite particular since, though it 
is the lowest moment, which a priori can be more sensitive to the lowest meson mass, its stability is pushed to 
higher energies because of the relatively large value of tpsiO)- As a biproduct, ci ^'^ becomes less sensitive to a 
low-mass resonance than . 



4.3 Mass of the gluonium G 

In order to study the properties of the gluonium G expected from other approaches to be in the range of 1.4 to 
2 GeV, we have to work with sum rules that are more sensitive to the high-energy region than . Using the 
positivity of the spectral functions, an upper bound on the gluonium mass squared can be obtained from the 
minimum (or inflexion point) of the ratios R: 

7^2l ^ -fjy, (27) 

giving (see Fig. 2): 

Mg < (2.16 ±0.16 ±0.14 ±0.05) GeV, (28) 

where the errors come respectively from the value of A, from the estimated unknown higher-order terms and 
from the gluon condensate. At such a small value of r, where the sum rule is optimized, we expect that high- 
dimension terms including instanton effects are highly suppressed. Combining these errors in quadrature, we 
obtain: 

Mg < (2.16 ±0.22) GeV, (29) 



Mq [GeV] 




Figure 2: r behaviour of the ratio of moments T^-jx . Its minimum corresponds to an upper bound for the 
scalar gluonium mass. 



An estimate of the mass squared can be obtained using the duality ansatz parametrization of the spectral 
functions, which leads to the FESR-like ratios 

dt e exp(-tT) i Im^G(i) , ^ ^ 

^21 = —f A^G, 30 

dt t exp(-tr) i ImV'G(i) 

where tc is the QCD continuum threshold. Neglecting in the analysis the eventual contribution of the light ub 
of a mass below 1 GeV, which is justified by the higher power of mass suppression in the present sum rule, we 
show in Fig. 3a, b the stabilities in t (here it is an inflexion point) and in tc (here it is a minimum) of the 
estimate. 



We have also studied the ratio 7?.io, but found no t stabihty there. 
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Figure 3: a): r behaviour of the estimated value of the scalar gluonium mass for — 4.5 GeV'^. The optimal 
value is at the inflexion point corresponding to the zero of its second derivative in r. b): behaviour of Mq. 

Then, we deduce the optimal result: 

Mg (1.50 ± 0.10 ± 0.15 ± 0.05) GeV, (31) 

where the errors are due to A, the estimated unknown higher-order terms and the gluon condensate. Combining 
these errors in quadrature, we obtain: 

Mg ^ (1.50 ± 0.19) GeV, (32) 

in agreement with the most recent sum rule result in |^ . Within the sum rule approach, one can approximately 
identify the value of tc as the mass squared of the next radial excitation: 

Mg' ^yTc^ (2.0 - 2.1) GeV, (33) 

where it can be noticed that, unlike the usual hadrons, the splitting between the lowest ground states and the 
radial excitations are relatively small {AMg/Mg ~ 30% compared with AMp/Mp > 70%). One can compare 
our value of Mg with the theoretical estimates (lattice calculations |^"||ll|, QCD inequalities |l2j, ...) and with 
the GAMS § and Crystal Barrel [|l| data. 

4.4 Decay constant of the G(1.5) 

The decay constant of the G(1.5) can be introduced via: 

me^^lG) = V2fGM^, (34) 

(1 2) 

where /g is the decay constant analogue of /jr. One can either estimate /g, by using the previous Cs ' sum 
rules, or using the SSR ci~'^\ as in NV, which we shall discuss later on. Using the previous sum rules, for 
instance ci^\ we obtain (see Fig. 4): 

/g (390 ± 98 ± 65 ± 39 ± 76) MeV, (35) 

where the errors come from A, the estimated unknown higher-order terms, the gluon condensate and the value 
of tc- Combining these errors in quadrature, we obtain: 

/g =i (390 ± 145) MeV, (36) 

which we shall use later on. However, one should remark that the errors in the present determinations are 
relatively large compared with the typical 10-20% accuracy in the QSSR estimate of the corresponding quantity 
in the ordinary hadron sector. 
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4.5 Mass and decay constant of the three-gluon bound state 

Recall that a QSSR analysis of the two-point correlator associated to the scalar three-gluon local current to 
leading order @, |l§: 

ha = g'fabcG^G'G' (37) 

leads to the mass prediction: 

Msg ^ 3.1 GeV y/Tc ~ 3.4 GeV, (38) 
and to the value of the decay constant of: 

/3G ^ 62 MeV, (39) 

which is relatively high compared with the mass of the gluonium built from the two-gluon current, and which 
makes the mass- mixing between these two gluonia states tiny ^: 

03 ~ 4° , (40) 

where the accuracy for the above predictions are about 10-20%. This state might be produced in the radiative 
decay of heavy quarkonia, when phase space permits, while its experimental search can be done by measuring 
some typical gluonia decays into the pairs ri'i]' or 7777'. 

4.6 Decay constants of the as and a'^ mesons 

As discussed previously and in NV, one can expect that the low moments ci ^'^"^ are sensitive to the low-mass 

(12) 

resonances whose effects can have been missed in the previous analysis of CI , and presumably in the lattice 
calculations within a one-resonance parametrization. In the following, we shall therefore test the "gluonium" 
nature of the broad low-mass states ctb and a'g, where the former, which has a mass in the range 0.5 to 1 GeV, 
is the one seen in 77 and tttt scattering experiments, and expected from the linear a model |4|], while we 
identify the second as its radial excitation with a mass close to that of the observed state at 1.37 GeV. We use 
a three-resonance {ctb, o'b G(1.5)) parametrization of the spectral function 0. We introduce the previous 
values of the G(1.5) parameters and the corresponding value of tc — 4.5 GeV^. Using as input Ma-'^ w 1.37 GeV 
and Mag ~ (0.5 ~ 1) GeV, we extract from the two sum rules the decay constants of the as and a'g within the 
r-stability criteria. We show in Fig. 5 the r behaviour of: 

h ^ fas (1 + P exp (M^^ - M^^ )r) , (41) 
h^faJl + Pj^ exp (Mj^ - M2 Jr 1 , (42) 



and of: 



where: 



<7b J '^B 



^Unfortunately, we cannot fix simultaneously the masses and decay constants of the erg and a'„ with the two sum rules. 
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|Mga= 1 GeV] 

U [GeV] 




Figure 5: t behaviour of the decay constants /2 from Cs and fi from Cs 



Then, we deduce: 



1.0 
1.4 
1.9 



/a' [GeV] 



0.6 
0.7 
0.5 



for AU^ [GeV] 



1.0 
0.75 
0.5 , 



(44) 



which indicates the necessity to have "low-mass gluonia" (which, as we shall see later on, couples strongly to 
tttt) for the consistency of the sum rules approach 0, although their effects are negligible in the high-moments 
analysis. The previous quantities will be useful later on for studying the as and a'g decays. 

5 Decay widths of the ctb-, ct'b Sind G 
5.1 as and a'^ couplings to tttt 



At this point, we use vertex sum rules to obtain further constraints. We consider the vertex (Fig. 6a): 
where: 



P1-P2 



1/(0) ^ 2mi 



(45) 
(46) 





11.11 



11.11 



b) 



Figure 6: Vertex controlling the giuonium couplings to: a) tttt; b) 77(77') — 77(77'). 



In the chiral limit {rn^ ~ 0), the vertex obeys the dispersion relation: 



dt 1, ^ 

-2 -luvV{t), 



which gives, by saturating with the three resonances crs, a'g and G: 

J X! 9S-^7rV2fs ~ 



(47) 



(48) 



S—(7b ,cr'p( ,G 



'^l^The existence of these states is also needed in the analysis of the correlator, by combining the LET result for ipsiO) with 
dispersive techniques [ po| . 
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Using the fact that V'{0) = 1 ^ from a generahzation of the Goldberger-Treiman relation using soft-pion 
techniques, one obtains a second sum rule: 

\ J2 9S..V2fs/Ml = \. (49) 

i=aB ,<y'g ,G 

Identifying the G with the G(1.5 ^ 1.6) at GAMS and the Crystal Barrel, we can neglect then its coupling to 
TTTT. Therefore, we can deduce: 

4 1 



V2/< 



173 



Fixing for definiteness M^r'^ ~ 1.37 GeV, one can deduce the width into tttt (tt+tt^, 27r°): 

r(cri3 ^ 7r7r)[GeV] w 0.7 r(crB ^ tttt) [GeV] w 1.3 for M^^ [GeV] w 1.0 

0.2 0.5 0.75 

0.1 0.9 0.5 , (51) 



where 



T{aB ^ tttt) = ( 1 _ ^] ^ . (52) 



We have repeated the derivation of /o-g by taking into account finite- width corrections. This leads to an increase 
of fas J which is compensated by the propagator effect in the estimate of gaBTnr, i-e. the result obtained from 
the vertex sum rule remains almost unchanged. It is interesting to see from the sum rule that a very light 
as around 500 MeV cannot be broad, which seems not to be favoured by the present data. Improvements of 
the data analysis are required for refining the mass measurement of the ctb. Our result indicates the presence 
of gluons inside the wave functions of the broad resonance below 1 GeV and the cr' (1.37), which can decay 
copiously into tttt |^. 

5.2 (7(1.5) coupling to r]r]' 

In order to compute these couplings, we consider the three-point function (Fig. 6b): 

%AQuq2) = fd*x, d^x2 e»fe-i+9--^) {0\TQ{xi)Q{x2)e^,{0)\0) , (53) 



where 9f_iu is the energy momentum tensor of QCD with three light quarks, while Q{x) is the topological charge 
density defined previously. Using the low-energy theorem: 

{m\o^^\m)^2Ml, (54) 

where rji is the unmixed U{1) singlet state of mass Alri^ — 0.76 GeV [Q, and writing the dispersion relation 
for the vertex, one obtains the NV sum rule: 

i ^ gs,nm^fs = 2Ml, (55) 

which implies, for Mo-^ ~ (0.75 ~ 1) GeV, and by assuming a G-dominance of the vertex sum rule: 

gamm ~ (1-2 - 1-7) GeV. (56) 
Introducing the "physical" rj' and rj through: 

1]' ~ cosBprji — sill 9 pr/s 

1] ^ ainOprji + cosOpTjs, (57) 



^The decays of the physically observed states will be discussed later on. 
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where |5^, 54 

6lp~ -(18 ±2)° (58) 

is the pseudoscalar mixing angle, one obtains: 

r(G ^ W) =i (5 - 10) MeV. (59) 

The previous scheme is also known to predict (see NV and 1^): 

r = rc^^/rcr,^. ~ 0.22, ga^^^sindpgaw =^ Tg,,,, ^ (1-1 - 2.2) MeV. (60) 

compared with the GAMS data r ~ 0.34 ± 0.13. This result, can then suggest that the G(1.6) seen by the 
GAMS group is a pure gluonium, which is not the case of the particle seen by Crystal Barrel [Q. 

5.3 cr^(1.37) and (7(1.5) couplings to Air through aBCB 

Within our scheme, we expect that the At: are mainly initiated from the decay of pairs of as- A LET analogous 
to that in previous cases, can be written for the estimate of gcdBCB- Using 

{aB\e>;,\aB) = 2Ml^ (61) 

and writing the dispersion relation for the vertex, one obtains the sum rule: 

gsas.n^fs = 2Ml^. (62) 



- E 



We identify the a'g with the observed /o(1.37). Using the observed value of the width into 47r ^ and 
extracting the S'-wave part, which we assume to be initiated from 2(7 we deduce |^: 

T{a'B ^ (47r)5) ^ (46 ~ 316) MeV =^ 9^'^.^^^ - (2-3 ^ 5.9) GeV, (63) 

where we have taken finite width corrections from the Breit-Wigner parametrization of the as resonance. We 
have used Mo-^ ~ 1 GeV and T{aB — > tttt) ~ 0.7 GeV. Neglecting, to a first approximation, the ctb contribution 
to the sum rule, we can deduce: 

gaasau ~ (2.7 ~ 4.3) GeV, (64) 

which, despite the large error, is larger than the Grjirji coupling. It can indicate that the decay of a pure 
gluonium state can be dominated (phase space permitting) by the An branching ratio initiated from the pair of 
cr-mesons, as already emphasized by NV. Using the previous values and the corresponding decay constant and 
width, one can deduce, using a Breit-Wigner form: 

r(G ^ ffBCTB An) ~ (60 - 138) MeV. (65) 

This feature seems to be satisfied by the states seen by GAMS and the Crystal Barrel. Our approach shows 
the consistency in interpreting the G(1.6) seen at GAMS as an "almost" pure gluonium state (ratio of the rj-q' 
versus the 7777 widths) , while the state seen by the Crystal Barrel, though having a gluon component in its wave 
function, cannot be a pure gluonium because of its prominent decays into 777/ and tt+tt". We shall see later on 
that the Crystal Barrel state can be better explained from a quarkonium-gluonium mixing. 

5.4 (Tb, cy's cind G couphngs to 77 

The two-photon widths of the cr, a' and G can be obtained by identifying the Euler-Heisenberg effective 
Lagrangian (Fig. 7a) | |5l| : 

C^q = —-rz 77[28Fui'Fi,\G\aGan + lAFnjyGi^xFxaGaa 

180?7i^ L ^ 

— Ffi^G ^vFapGafj — F^ijF^^GapGap], (66) 



^^We have taken the largest range deduced from the different branching ratios given by PDG [Q. We have used T{cfb — ► all ~ 
(150 ~ 500) MeV. 
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where Qq is the quark charge in units of e, —f3i =9/2 for three flavours, and is the "constituent" quark 
mass, which we shall take to be 



with the scalar-77 Lagrangian 



niu ~md ~ Mp/2, nis ~ M$/2 , 



(67) 
(68) 





(a) (b) 
Figure 7: Vertex controlling the gluonium couplings to J/^p(-^)^: a) box diagram; b) anomaly diagram. 



This leads to the sum rule: 

q=u,d^s 

from which wc deduce the couplings 

g^^^^ « (0.4 - 0.7)a GeY'^ g„'^^^ « (0.5 - 0.6)a GeV"^ 

Using the corresponding decay width: 



T{S 77) = 



|gS77l' 

167r 



M: 



Si 



one obtains: 

r(crB ^ 77) « (0.2 - 0.3) keV r(crB ^ 77) « (0.7 - 1.0) kcV 



(69) 

ffG77 ~ (0.6±0.2)a GeV"^ (70) 

(71) 

r(G ^ 77) « (1.0±0.8) keV, (72) 



in agreement with the NV results, but smaller (as expected from general grounds) than the well-known quarkonia 
widths: 

r(^' 77) - 4.2 keV r(/2 ^ 77) ~ 2.6 keV. (73) 
Alternatively, one can use the trace anomaly (Fig. 7b): 



(0|^;:i7i72) = (0|i/3(«.)G^ + gFri^ri7i72), 



(74) 



where F^^^ is the photon field strength and R = Sj^Qi- Using the fact that the RHS is 0{k'^) implies the sum 
rule g 0: 



-/3(a,)G2|7i72) 



. aR 



FrFr hii2), 



from which one can deduce the coupling: 



V2 



X! fs9s 



S — (TB -X^'n -G 



aR 
3^- 



(75) 



(76) 



It is easy to check that the previous values of the couplings also satisfy the trace anomaly sum rule. However, 
the result for the (777 coupling is much smaller than the ones in [ p7| , ^ obtained from a single resonance 
saturation of the trace anomaly, because (as we have seen previously) the one resonance saturation is not a 
good approximation, while the value of f^^^ used in these papers is much smaller than the one obtained here. 



*Here and in the following, wc shall use AI^ 



(0.75 - 1.0) GeV. 
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5.5 cxb, ct'b and G productions from radiative J/ip decays 

As in | ]5l[ , one can estimate this process, using dipersion relation techniques, by saturating the spectral function 
by the J/^p plus a continuum. The glue part of the amplitude can be converted into a physical non-perturbative 
matrix element (0|asG'^|S') known through the decay constant fs estimated from QSSR. By assuming that the 
continuum is small, one obtains: 



/3?656100 V Mc J \Mc) r(J/V' ^ e+e") 

where Mc — 1.5 GeV is the charm constituent quark mass. We use — /3i = 7/2 for six flavours. This leads to 
the rough estimates: 

.B(J/?/> -> 7(75) X B(crB ^ all) « (0.4- 0.6) X 10-3 
B(J/^ -> 7(7^) X ^((7^ ^ all) w (0.8 - 1.0) X 10^3 
B( J/V' ^ iG) X B{G all) w (0.5 - 0.4) x 10"^ (78) 

These branching ratios can be compared with the observed B{J/ip jrj') and B{J/tfj if 2) ones, which are 
respectively 4 x 10-^ and 1.6 x 10"'^. The ctb could already have been produced, but might have been confused 
with the TTTT background. The "pure gluonium" G production rate is relatively small, contrary to the naive 
expectation for a glueball production. In our approach, this is due to the relatively small value of its decay 
constant, which controls the non-perturbative dynamics. Its observation from this process should wait for the 
rCF machine. However, we do not exclude the possibility that a state resulting from a quarkonium-gluonium 
mixing may be produced at higher rates. From the previous results, one can also deduce the corresponding 
stickiness defined in 

6 Properties of the scalar quarkonia 

6.1 Mass and decay constant of the S2{uu + dd) quarkonium 

We consider this state as the SU{2) partner of the ao(0.98) associated to the divergence of the charged vector 
current of current algebra: 

d.yf'ix) = (m„ - md)uiij5)d. (79) 
The mass and coupling of the uq have been studied within the QSSR [l^ : 

Ma, ~ (1 - 1.05) GeV fao ^ (0.5 - 1.8) MeV, (80) 

where the small value of fao is due to the light current quark mass difference. In our approach, due to the good 
realization of the SU(2) symmetry, the mass of the 6*2 (uu + dd) bound state is expected to be degenerate with 
the one of the ao(u(i). The continuum threshold at which the previous parameters have been optimized can 
roughly indicate the mass of the next radial excitation, which is [|l6| : 

Ms- « V^T - (1.1 - 1-4) GeV » M^,, (81) 

which is about the /o(1.37) mass. An estimate of this mass using a model with an infinite number of resonances 
leads also to the same result (see e.g. S.G. Gorishny et al. in [^). 

6.2 Couplings of the 5*2 to tt+tt , K^K^ and 77 

Using vertex sum rules, BN | |2^ ] obtain the ^2 coupling to pair of pions in the chiral limit: 

gs,^^^- ^ ^ (uu)re*^^ * 2.46 GeV, (82) 

for the typical value of r = 1 GeV^^, in good agreement with the one from the SU{?>) relation between the 
5*2 — > TTTT and oq rjn widths and with gaoK+K- — 2 GeV from [^ , as intuitively expected. Therefore, we 
deduce with the same accuracy as the one of the measured qq — > rjn width: 

r(S'2 ^ TT+TT") ~ 120 MeV. (83) 
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Using SU{3) symmetry, one can also expect: 



1 

gs2K+K- - -^gs2-^+n- ■ (84) 

The 77 width can also be obtained from the prediction in BN: 

25 

r{S2 -> 77) - yr(ao(.98) -> 77) ~ 0.67 keV, (85) 

where 25/9 is the ratio of the u and d quark charges, the last number coming from the data. This result is 
supported by the vertex sum rule analysis ||5^ . 

The estimate of the 77 and hadronic widths of the 52 is more uncertain. Using the phenomenological observation 
that the coupling of the radial excitation increases as the ratio of the decay constants r = fs^/fs'^i we expect: 

3 



r(5^-77) « ns2^ii), 



V{S'^^nn) « r^i^y^S^^^^), (86) 

which by taking r (M^^/Msj)*^"^^^, like in the pion case gives |^ 

77) « (4 ± 2) keV, T{S2 TT+TT") w (300 ± 150) MeV. (87) 

To a first approximation, we expect that the decay of the S2 into At: comes mainly from the pair of p mesons, 
while the one from 2aB (gluonia) is relatively suppressed as aj. using perturbative QCD arguments. 

6.3 Mass and decay constant of the 6*3(55) quarkonium 

In order to complete our discussions in the scalar sector, we compute from the sum rules the mass and decay 
constant of the S^{ss) state. In so doing, we work with the two-point correlator: 

i;-Uq^) = 1 j d^x e*«^ {Q\T.h,{x) (J,.(0))t |0), 

where: 



Jss{x) = nisss, (89) 



and we introduce the S3 as 



{0\J,s\S3) ^ ^fssMl. (90) 
We work with the Laplace transform sum rules: 



f°° 1 
^ss{t) = / di exp (-tr) -Im7/'3(i), 

J 4m? I" 



/ 4m? 



^ _ Ximsdi i exp(-tT) ilmi/)j^(i) 

n., M, 



2 



-'"if* (1.43) 



(91) 



The QCD expressions of the sum rule have been obtained in [ ]60| and are now known to three-loop accuracy 
(see the compilation in |6^). A much better stability is obtained at r ~ 0.4 ~ 0.6 GeV^ by working with 
the double ratio of sum rules instead of the ratio. Using mi (lGeV) ~ (150 ^ 190) MeV correlated to the values 
of A [Q, we deduce (see Fig. 8a, b): 

MgjAfj^.(i.43) ~ 1.03 ± 0.02 =^ Mgs ~ (1474 ± 44) MeV. (92) 



^^We estimate the error by assuming that n = 2 ± 1 
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This result confirms the earher QSSR estimate in @ 0- 




Figure 8: r and tc behaviours of M5s/M;f»(i 43). 



The resuh indicates the mass hierarchy: 

Ms, ^uu+M < Mk; ^us < M,s . (93) 

The SU{S) breaking obtained here is shghtly larger than the nai've expectation as, in addition to the strange- 
quark mass effect, the (ss) condensate also plays an important role in the splitting. The sum rule J-gg leads to 
the value of the decay constant (Fig. 9): 

= (43 ± 19) MeV. (94) 




These previous estimates have been optimized at tc — 6 GeV^. Therefore, we expect that the radial excitation 
S'g will be in the range: 

Ms^ « (1.7 ~ 2.4) GeV, (95) 

where the first number corresponds to the phenomenological extrapolation AIs;^ — Ms, — Ms' — Ms^, while the 
second value is -v/t 



6.4 Couplings of the S3 to K^K and 77 

In so doing, we work with the vertex (Fig. 10): 



V{p,q) = / d'x d'y e-^f-e*(''+^)^(0|TJ^+(x)JK-(y)J,.(0)|0) , (96) 



and evaluate it at the symmetric point = = [p + qY ~ ~iQ^ > 0) as in p8{ . 

^® After obtaining this result, we learned that a similar value has been obtained from lattice calculations | ]m| . 
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+ . . . 



q q 

Figure 10: Vertex sum rule for the Sj,K'^K~ coupling. 



Then, we obtain, in 4 — e dimensions: 



3 fQ'^\^^^^ 1 _ 

__,,,r(./2)(^— j +—{uu^ss) 



Its phenomenological expression can be approximated by: 



2/|m4 2flMi 



7T\9S3K+K- 



The Laplace transform of the previous equation leads to: 

9S,K+K- V2Mlfss ^ xdx e-(^?=-^^l:)- ~ 4ms - y t(wu - 2, 
Using [pl (ss)/(uu) ~ 0.7, (uu){l GeV) ~ -(229MeV)3 and m, (1 GeV) ~ (150-190) MeV, we obtain: 



(97) 



(98) 



(99) 



9SiK+K- - (2.7 ± 0.5) GeV and gsaw - 0.9W -Qs^k+k- ■ 



We can also predict: 



9S3 



sm Upgs^ 



(100) 



(101) 



This result can be compared with the one obtained previously, and with gagK+K- — 2 GeV from Ipq] , which 
expresses the good SU(3) symmetry of the couplings as intuitively expected. Therefore, we can deduce: 



T{S3 K+K-) ~ (73 ± 27) MeV and TiSs -> 7777) ~ (15 ± 6) MeV. 



(102) 



A comparison of this result with the /o(1.5) width into KK of 5 MeV, and the strong coupling of the /o(1.5) 
to TTTT leads us to the conclusion that this state cannot be a pure ss state. 

We estimate the 5377 width using its relation with the one of 6*2 and the corresponding quark charge: 



r(S'3 ^ 77) w 4 (tt^^ r(S'2 77) ~ (0.40 ± 0.04) kev. 



25 \Ms, 

Finally, analogous to the case of 5*2, we can also have for Mg/^ ~ 1.7 GeV: 

^{S's 77) ~ (1.1 ± 0.5) keV, r(S'3 K+R-) w (112 ± 50) MeV. 



(103) 



(104) 



7 "Mixing-ology" for the decay widths of scalar mesons 

We have discussed in the previous sections the masses and widths of the unmixed gluonia and quarkonia states. 
The small value of the mass mixing angle computed in Q from the off-diagonal two-point correlator (see Fig. 
11), which is proportional to the light quark mass, allows us to neglect the off-diagonal term in the mass matrix, 
and to identify the physical meson masses with the ones of the unmixed states. 
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Figure 11: Off-diagonal two-point correlator controlling the quarkonium-gluonium mass mixing angle. 



In the following, we shall be concerned with the mixing angle for the couplings, which, in the same approach, 
is controlled by the off-diagonal non-perturbative three-point function (see Fig. 12) which can (a priori) give a 
large mixing angle. 




+ 



Figure 12: Off-diagonal non-perturbative three-point function controlling the coupling mixing angle in the 
quarkonium-gluonium decays into pion pairs. 



7.1 Mixing below 1 GeV and the nature of the a and /o(0.98) 



This part will be an update of the scheme proposed by BN |27 . We consider that the physically observed /o 
and (T states result from the two-component mixing of the as and 5*2 = "^(^'^ + '^d) unmixed bare states: 



l/o) = -sines\<JB)+cos9s\S2) 
\a) = cos 9 s\(J b) +sva0s\S2) 



(105) 



We also use the previous prediction for T{aB — * 77) — (0.2 ~ 0.3) keV, and the experimental width r(/o 
77) « 0.3 keV. Therefore, we can fix the mixing angle to be: 



(40- 45)°, 



(106) 



which indicates that the /o and cr have a large amount of gluons in their wave function. This situation is 
quite similar to the case of the t]' in the pseudoscalar channel (mass given by its gluon component, but strong 
coupling to quarkonia). Using the previous value of 6*5, the predicted value of gs^K+K^-, the approximate 
relation gs^K+K- — ^feir+u--' ^'^'^ the almost universal coupling of the gb to pairs of Goldstone bosons, one 
can deduce: 



5/o7r+7r- - (0.1 2.6) GeV, gfoK+K- 

garr+rr- - g^K+K- ^ (4 5) GcV, 



-(1.3 - 4.1) GeV, 



(107) 



which provides a simple explanation of the exceptional property of the fo (strong coupling to KK as observed 
in TTTT and KK data ^^), without appealing to the more exotic four-quark states and KK molecules |^. 
Using the previous predictions for the couplings, and for 63, we obtain for Afg-B — (0.75 — 1.0) GeV: 



r(cr ^ TT+TT-) « (300 - 700) MeV, 



(108) 



r(/o(0.98) ^ TT+TT") « (0.2 ~ 134) MeV 

^^A OSSR, analysis of the ao(0.98) within a four-quark scheme leads to too low a value of its 77 width as compared with the 
data [BSl. 
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By recapitulating, our scheme suggests that around 1 GeV, there are two mesons that have 1/2 ghion and 1/2 
quark in their wave functions resulting from a maximal destructive mixing between a quarkonium (52) and 
gluonium (ct_b) states: 

• The (ao(0.98)) is narrow, with a width < 134 MeV, and couples strongly to KK, with the strength 

9 foK+ K~ / 9 fo-n+TT- — 2- This property has been seen in tttt and 77 scatterings and in pp experiments, 
and also suggests its production from the radiative decay. 

• The CT, with a mass around 0.75 ^ 1 GeV, is large, with a width of about (300 ~ 700) MeV. 

7.2 Mixing above 1 GeV and the nature of the /o(1.37), /o(1.5) and /j(1.7) 
7.2.1 The data 

Let us recall the experimental facts |^. Coupled channels analysis demonstrate that in addition to the broad 
cr(l.O) below 1.2 GeV, one needs 0the /o(1.37), with the following properties 

r(/o(1.37) all) ~ (150 - 500) MeV r(/o(1.37) ^ 77) ~ (5.4 ± 2.3) keV, (109) 

(where the number for the 77 width is less reliable); however, from the quoted branching ratios and extracting 
the All S-wave, we deduce psf : 

r(/o(1.37) ^ tttt) < (15 - 100) MeV, r(/o(1.37) ^ (47r)5) ~ (46 - 316) MeV. (110) 

For the /o(1.5), one has @ |: 

r(/o(1.50)-. all)c.(100^150)MeV, ^ (3.4 ± 0.8), (111) 

r(/o(1.50) ^ tttt) 

while the partial widths of the /o(1.5) divided by the phase-space factor satisfy the ratios: 

TTTT : 777; : W : KK~1 : (0.27 ±0.11) : (0.19 ±0.08) : (0.24 ±0.09). (112) 

By assuming that the branching ratios other than the previous ones are negligible, we can deduce the experi- 
mental data: 

r(/o(1.50) ^ tttt) ~ (20 - 31) MeV r(/o(1.50) ^ KK) ~ (3.6 - 5.6) MeV 

r(/o(1.50) -> r/?7) ~ (2.6 -3.3) MeV r(/o(1.50) -> 47r°) ~ (68 - 105) MeV, (113) 

and: 

r(/o(1.50) W) - 1-3 MeV. (114) 

In the region between 1 and 1.5 GeV, we shall be concerned with the four unmixed states: the radial excitations 
S'^(1.3) of the 52 and crs(1.37) of the as, the S'3(1.47) and the G(1.5). 



7.2.2 Nature of the /o(1.37) 

By inspecting the experimental data in the region around 1.37 GeV, one can expect that the /o(1.37) which 
has a large 77 width, is a good candidate for being essentially composed by the quarkonium 5*2. However, 
its large decay width into An indicates its large gluon component from the a'^. Due to the almost degenerate 
values of the two previous unmixed states, we expect that they will mix in a maximal way, as in the case of 
their corresponding ground states 5*2 — ctb of mass around 1 GeV. Therefore, we expect that it is difHcult to 
disentangle the effects of these two resonances where: 

• The S2 quarkonium state is (relatively) narrower with a total width around 400 MeV. It decays mainly into 
TTTT (its decays into KK and 7777 can be obtained by taking proper Clebsch-Gordan factor and assuming SU{3) 
symmetry of the couplings.), and has a large 77 width of about 4 keV. 

• The a' is very broad, with a width as large as 1 GeV, has an universal coupling to tttt, KK and 7777, and can 
decay into An via aa, with a width of (46 ~ 316) MeV. However, its decay into 77 is relatively small. 

• The experimental candidate /o(1.37) has amusingly the combined properties of the 6*2 and a' . 



In the following we shall take the largest range of values deduced from the ones (sometimes controversial) given by PDG 
The question whether there is one or two /o(1.37) (one coupled strongly to 2it and the other to 4it) is not yet settled. 
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7.2.3 A 3x3 mixing scheme and nature of the /o(1.5) 

In the following estimate of the hadronic widths of the remaining other mesons, we expect that, due to its large 
width, the a' can have a stronger mixing with the 5*3(1.47) and the G'(1.5), than the 5*2, such that, to a first 
approximation, we shall only consider the mixing of the three states (t'(1.37), 5*3(1.47) and the G(1.5) through 
the unitary 3x3 CKM-like mixing matrix of the couplings 



/o(1.37) \ / C12C13 S12C13 si3 \ / (t'(1.37) 

/o(1.50) = -S12C23 - C12S23S13 C12C23 - S12S23S13 S23C13 53(1.47) I (115) 

/o(1.60) / \ S12S23 - C12C23S13 -C12S23 - S12C23S13 C23C13 / \ 0(1.5) 

where: 

Cij = cos 9ij and = sin^^. (116) 
We shall use in our numerical analysis: 

r(cr' tttt) w (1-^2) GeV T{ct' ^ in) w (40 - 316) MeV (data), 

T{(7'^KK) w (0.7 -1.5) GeV r(cr' r/ry) w (0.2 - 0.4) GeV, (117) 

and we shall discuss separately the cases of large (upper values) and small (lower values) of each partial widths. 

• The case of small widths 

We estimate each entry as follows: 

For the first line of the matrix, we shall fix C12 from the small width of /o(1.37) into KK, which gives: 

ci2~0.31 si2«-0.95. (118) 
Then , we fix ^13, by requiring the best prediction: 

r(/(f (1.37) ^ 47r) w 150 MeV (119) 
compared with the previous data. This implies the two solutions (a) and (b): 

ci3 ~ (a) - 0.29 and (b) 0.74. (120) 
For the second line of the matrix, we use the observed width r(/o(1.5) — > tttt), in order to get: 

C23 ~ (a) 0.45 and (b) 0.37. (121) 

for the corresponding two values of C13. Finally, the observed width r(/o(1.5) — > Att) favours the alone case (b). 

• The case of large widths 

We repeat the previous analysis taking the upper values of the partial widths. 

• Predictions 

Our final mixing matrix is given by the largest range spanned by each mixing angles. Therefore, the mixing 
matrix reads: 

/o(1.37) \ / 0.01 - 0.22 -(0.44- 0.7) 0.89 - 0.67 \ / cr'(1.37) \ 
/o(1.50) w 0.11 - 0.16 0.89 - 0.71 0.43- 0.69 53(1.47) , (122) 
/o(1.60) / \ -(0.99 - 0.96) -(0.47 - 0.52) 0.14- 0.27 / \ G(1.5) / 

where the first (resp. second) numbers correspond to the case of large (resp. small) widths. From the previous 
schemes, we deduce the predictions 0: 

r(/o(1.37) ^ tttt) w (22 - 48) MeV r(/o(1.37) ^ w) < 1. MeV r(/o(1.37) ^ W) < 2.5 MeV (123) 



^"in the following, we shall neglect the CP-violating phase. It can also be noticed that we cannot fix the masses of the physical 
states from this mixing of the couplings. However, as shown earlier from the small mass-mixing angle, we expect that their masses 
are around those of the unmixed states. 

^^Recall that we have used as inputs: r(/o(1.37) — > KK) ft: 0, r(/o(1.5) — > tttt) fa 25 MeV, while our best prediction for 
r(/o(1.5) — > (47r)g) is about 150 MeV. The present data also favour negative values of the foW^ fov'v ^^'^ foKK couplings. 
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and 



r(/o(1.5) KK) « (3 ~ 12) MeV r(/o(1.5) ^ 7777) « (1 ~ 2) MeV r(/o(1.5) ^ W) < 1- MeV . (124) 

Despite the crude approximation used and the inaccuracy of the predictions, these results are in good agreement 
with the data, especially from the Crystal Barrel collaboration. These results suggest that the observed /o(1.37) 
and /o(1.5) come from a maximal mixing between the gluonia (cr' and G) and the quarkonium S'3 states. The 
mixing of the S'3 and G with the quarkonium 6*2, which we have neglected compared with the u', can restore 
the small discrepancy with the data. One should notice, as already mentioned, that the state seen by GAMS is 
more likely to be similar to the unmixed giuonium state G (dominance of the 47r and 7777' decays, as emphasized 
earlier in NV), which can be due to some specific features of the central production (double pomeron exchange 
mechanism which favours the gluonia production) for the GAMS experiment. This feature is not shared by the 
data from the Crystal Barrel and Obelix collaborations, which correspond to the annihilation of antiprotons at 
rest in a hydrogen target, and therefore can also favour the production of quarkonia states |^ 



7.2.4 Nature of the /./(l. 71) 
For the /o(1.6), we obtain: 

r(/o(1.6) -> KK) K (0.5 - 1.6) GeV r(/o(1.6) ^ tttt) « (0.9 - 2.) GeV 

r(/o(1.6) ^ 7777) « (0.04 - 0.6) GeV r(/o(1.6) ^ t^t;') » (0.03 - 0.07) GeV, (125) 

and 

r(/o(1.6) -> (471)5) « (0.02 ~ 0.2) GeV, (126) 

which suggest that the /o(1.6) is very broad and can again be confused with the continuum. Therefore, the 
/j(1.7) observed to decay into KK with a width of the order (100 ~ 180) MeV, can be essentially composed 
by the radial excitation 5*3(1.7 ~ 2.4) GeV of the S^{ss), as they have about the same width into KK (see 
section 6.4). This can also explain the smallness of the /j(1.7) width into tttt and 47r. Our predictions of the 
/j(1.71) width can agree with the result of the Obelix collaboration while its small decay width into Air is 
in agreement with the best fit of the Crystal Barrel collaboration (see Abele et al. in |lj]), which is consistent 
with the fact that the /o(1.37) likes to decay into Air. However, the /o(1.6) and the /j(1.71) can presumably 
interfere destructively for giving the dip around 1.5 ^ 1.6 GeV seen in the KK mass distribution from the 
Crystal Barrel and pp annihilations at rest []68| , |67| . 



7.2.5 Comparison with other scenarios 



One can also compare our results with some other mixing scenarios [ |69| , 70 1. Though the relative amount of 
glue for the /o(1.37) and /o(l-5) is about the same here and in ^ , one should notice that, in our case, the tttt 
partial width of these mesons come mainly from the cr', a glue state coupled strongly to the quark degrees of 
freedom, like the rj' of the ^7(1)^ anomaly, while in |69j, the S2 which has a mass higher than the one obtained 
here plays an essential role in the mixing. Moreover, the /j(1.71) differs significantly in the two approaches, as 
here, the /j(1.71) is mainly the ss state S'3, while in []69| , it has a significant gluon component. In the present 
approach, the eventual presence of a large gluon component into the /j(1.71) wave function can only come from 
the mixing with the broad /o(1.6) and with the radial excitation of the giuonium G(1.5), which mass is expected 
to be around 2 GeV (value of ^/t^ in the QSSR analysis). However, the apparent absence of the /j(1.71) decay 
into 47r from Crystal Barrel data may not favour such a scenario. 



7.2.6 Summary 

Within the present approach and the present data, we expect above 1 GeV that: 

• the /o(1.37) is a superposition of two states: the radial excitation S2 of the quarkonium S2 = {uu + dd) and 
an /o(1.37) resulting from a maximal mixing between the giuonium G and the S^^ss) with the radial excitation 
a'g of the broad low mass <7b- Its large 77 width comes from the S2, while its afhnity to decaying into An comes 
essentially from the cr'. 

• the /o(1.5) with the properties observed by the Crystal Barrel collaboration results from a maximal mixing 
between the giuonium G and the S3(ss) with the radial excitation cr' of the broad low-mass cr. This large gluon 

^^We plan to come back in more details to this point in a future work. 
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component explains its affinity to decaying into 27r (from (t'(1.37)) and An (from a' and G) 
• the /j(1.71) can be identified with the radial excitation S'^ of the ground state 53(1.47), as they have about the 
same width into KK. This can also explain the absence of its decay into tttt and 47r. The eventual observation 
of these decays can measure its expected tiny mixing with the wide (and presumably unobservable) /o(1.6) and 
with the radial excitation of the G{1.5). The presence of the dip around (1.5 ~ 1.6) GeV seen in the KK 
invariant mass distribution may already be a signal of a such mixing. However, we need a clear spin-parity 
assignement of the /j(1.71) before a definite conclusion on its nature can be drawn. 



8 The tensor gluonium 
8.1 Mass and decay constant 

In this section, we shall study the properties of the 2++ gluonium. We shall be concerned with the two-point 
correlator: 

VtJ.p'Hu'y + Vt^oVup - T,lt^ur]pa] i'riq'^), (127) 



where: 



3 



Vt^y = g^iu (128) 



To leading order in as and including the non-perturbative condensates, the QCD expression of the correlator 
reads 0: 

i,T{q^ ^ -Q2) ^ _ 1 (g4) log % + l|^(20i - O2) , (129) 



where: 

Using the vacuum saturation hypothesis, one can write: 



0\ — ifabcGtiaGjya)^ and O2 — UabcG^uGap f' ■ (130) 



{20,-02)c,-^{G^f. (131) 



In order to get the gluonium mass, we work with the ratio of moments: 



£(1) 

Tl,o ^ (132) 



as in SN Its QCD expression reads: 



nw^3T-'\l + ^7r'^{asGyy (133) 



4.0 - 




Figure 13: r behaviour for the upper bound of the tensor gluonium mass AIt- 
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Using the positivity of the spectral function, the minimum of the ratio of the moments (Fig. 13) leads to the 
upper bound: 

Mt < (2.7 ± 0.4) GeV, (134) 
while a resonance + QCD continuum parametrization of the spectral fimction gives the value (Fig. 14): 

Mt (2.0 ± 0.05 ± 0.05 ± 0.05) GeV ~ (2.0 ± 0.1) GeV. (135) 




5 2-0 



1.9 



0.05 O.IO^aiS 0.20 0.25 0.30 
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4.5 5.0 5.5 6.0 
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Figure 14: a) r and b) tc behaviours of the 2++ tensor gluonium mass. 

The mass obtained here is larger than the one in SN, which is mainly due to the increase of the gluon condensate 
value. The errors come respectively from the gluon condensate, the factorization assumption and the continuum 
threshold tc- The decay constant can be extracted from and reads (Fig. 15): 

/t =i (80 ± 8 ± 4 ± 4 ± 10) MeV ~ (80 ± 14) MeV, (136) 

where the first error comes from Mt and the remaining ones are of the same origin as before. 
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Figure 15: r behaviour of the 2++ tensor gluonium decay constant /t. 

The value of tc at which these results are optimal is: 

y/Tc « Mt' (2.2 - 2.3) GeV, 



(137) 



which corresponds roughly to the mass of the radial excitation. Our result is slightly higher than the one in SN 
and in Q, because of the increase of the value of the gluon condensate. We do not expect that the perturbative 
radiative corrections will affect noticeably our prediction for the mass |^ , as this effect tends to cancel out in 
the ratio of moments. 



8.2 Tensor gluonium decay widths 

Using the result in one can constrain the ratio: 

gTTTTT 

r = 

5/2 

^^An investigation of these effects with K. Chetyrkin and A. Pivovarov is under way. A preliminary result indeed indicates that 
the effects are small. 



(0.68 ±0.23), 



(138) 
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using the data, where: 

nT^^-^-) = ^^\kA\ (139) 

being the pion momentum. Assuming an universal coupUng of the T to the pairs of Goldstone bosons, this 
leads to the width: 

r(r ^ TTTT + KK + ?7?7) < (119 ± 36) MeV, (140) 

which indicates that the 2++ cannot be wide, contrary to some claims in the literature. Another indication of 
the smallness of the 2++ width can be obtained from the low-energy theorem analogue of the one for the scalar 
gluonium used in section 5 ||71|. Using a dipersion relation, one can write: 



\ V2fT9T.. « 1. (141) 



4 

/2,T 

As injra, we assume that the vertex is saturated by the /2(1.27) and by the gluonium T. Using fj^ ~ 0.19 
GeV ]73[ |, the previous value of /t, and the experimental value gf^-jrv — 16 GeV~^, one can notice that the sum 
rule is already saturated by the /2-meson, and leads to: 

Sttttt ~ 2.6 GeV"^ r{T tttt) w 10 MeV, (142) 

while taking a typical 10% error in the estimate of ff^ , one can deduce: 

r(T ^ tttt) < 70 MeV. (143) 



The 77 width can be obtained by relating it to the one of the 0^+ gluonium within a non-relativistic quark 
model approach ^: 

r(T ^ 77) ^ ^ (^^^ r(0++ ^ 77) ^ 0.06 keV, (144) 
which shows again a small value typical of a gluonium state. 

8.3 Meson-gluonium mass mixing and the nature of the /j(1.71) and C(2-2) 

In order to evaluate the gluonium-quarkonium mass mixing angle, we shall work, as in 1^, ^ , with the off-diagonal 
two-point correlator: 



V^^V ^ rf^xe^«-(O|r0^,(x)0^,(O)t|O) 



where 

ei,{x) = iq{x){^^:D, + ^,D^)q{x). (146) 

Here, 13^ = 15^ — is the covariant derivative, and the other quantities have already been defined earlier. 
Taking into account the mixing of the currents, one obtains Q: 

^.(.^ - -Q^) - (^) ( V ^ - ?J log ^) _ log % ia.G% (147) 

IStt^ V tt / V 15 I 367r 



The resonance contribution to the spectral function is introduced using a two-component mixing formalism 
1 

TT 



-lmiljgq{t)~sin2eTM^fTMlfF5{t-M^) + "QCD continuum", (148) 



where 9t is the mixing angle; F is a generic notation for the quarkonia /2(1.27) and /2(1.53); M is the average 
of the G and F mass squared; Jt,f is the decay constant, where //j ~ (0.11 ~ 0.14)M/2 |7^. Noting that 
within our approximation, the Laplace transform sum rule does not present a r stability, we work with the 
FESR: 

£'.,l:„,„„.._|L(|)<3^,,|g(G»,), 
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where ^ lit'J + tf) ~ 2.2 GeV^ is the average of the continuum threshold for the quark and gluonia channels. 
Then, one obtains the approximate estimate |^ 0: 

Ot -10°, (150) 

which is a small mixing angle, and which indicates that: 

• the interpretation of the /j(1.71) as the lowest 2++ ground state gluonium, is not favoured by our present 
result. 

• the C(2-2) observed by the BES collaboration is a good gluonium candidate (mass and width), but it 
can be more probably the first radial excitation of the lowest-mass gluonium (tc ~ 2.2 GeV). Our result, for 
the mass, suggests that one needs further data around 2 GeV for finding the lowest mass 2++ gluonium. Its 
(non-)finding will be a test for the accuracy of our approach. Our prediction for the mass and for tc, which 
suggests a rich population of 2"'""'" states in this region, should stimulate a further independent experimental 
test of the gx candidates seen earlier by the BNL group in the OZI suppressed reaction n^p — + cjycjin [Q, and a 
further understanding of their non-observation in J/i/; radiative decays ||76|. 



9 Pseudoscalar gluonia 
9.1 Mass and decay constant 

The pseudoscalar gluonium sum rules have been discussed many times in the literature, in connection with 
the gluonium mass in SN 77' meson, the spin of the proton and the slope of the topological C/(l) charge 
in 1 77, Here, wc update the analysis of the gluonium mass in SN, taking into account the new, large 

perturbative radiative corrections [Q, and the radiative correction for the gluon condensate In so doing 
we work with the Laplace transform: 



L 



(0) = 



.00 

I dt exp(— ir) — Im?/'p(t) 

{a,G^)T^ (1 - exp(-t,T)) - 27rV3g(/,fe,G3,^) 



llTT 



(151) 



where: 



^pert 



P2 



1 + tcT 



Pi 

exp(-tcr) 



83 + 2/3i(l - 2-ie) - 4^ log(- logrA^) 

P2 



and the ratio of moments 



'(1) 



10 



d (0) Cp 

L,p 



(152) 



(153) 



Subtracting the 77' contribution^, and using the positivity of the spectral function, the minimum of 72,io gives 
the upper bound (Fig. 16): 

Mp < (2.34 ± 0.42) GeV, (154) 

where the error comes mainly from the value of A, which one can understand as the minimum occurring at 
large values of r. 

^*One should notice that, to this approximation, we do not have tc stabihty, such that our result should only bo taken as a crude 
indication, but not as a precise estimate. 



^^In order to take i 
constant obtained in 



account the change in the radiative correction coefficient, we have redone the estimate of the t]' decay 
[ we obtain a slight change f^/ 2; 30 MeV in agreement, within the error, with the previous value of 

2 _ ,,f2 _ 2 ,^2 _ (Q g7 GcV)2. 



(24 ± 3.5) MeV. We have used in the massless quark limit M^, = NP, 



3 if 
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Figure 16: r behaviour of the upper bound for the pseudoscalar gluonium mass Mp. 

The estimate of the pseudoscalar gluonia mass is given in Fig. 17, where the stabiHty in r is obtained at smaller 
T values, which gives smaller errors. We obtain: 

Mp ~ (2.05 ± 0.05 ± 0.13 ± 0.13) GeV ~ (2.05 ± 0.19) GeV, (155) 

where the first error comes from tc, the second and the third ones come from A and the truncation of the 
QCD series. As in the previous cases, we have estimated the unknown coefficient to be of the order of ±400 
(normalized to the lowest order term) by assuming a geometric growth of the QCD series. 




to[GeV2l 



Figure 17: r and tc behaviours of the pseudoscalar gluonium mass. 



The corresponding value of tc is: 

^/^c w Mp, ~ (2.1 ~ 2.3) GeV. (156) 

The decay constant has a good r stability (Fig. 18) though the tc stability is only reached at — 7 GeV^ (Fig- 
21). Using the range of tc values from 5 to 7 GeV^, we deduce: 

/p ~ (8 ~ 17) MeV. (157) 



fp [MeV] |tc = 5GeV2| 




tc[GeV2l 



Figure 18: r and tc behaviours of the decay constant fp for tc = 5.5 GeV^. 

These results are slightly higher than the one in SN obtained from a numerical least-squares fit, because mainly 
of the effect of the radiative corrections. 
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9.2 Testing the nature of the £'/i(1.44) 

Old and new experimental data indicate the presence of some extra states in the range of 1.40-1.45 GeV, 
which go beyond the usual nonet classification. In order, to test the nature of such states (denoted hereafter as 
the £'//,(1.44)), let us now come back to the sum rule: 

l'^p^^ = — exp(-tr) - IniV;p(t), (158) 

which has been used in for fixing the decay constant f,y : 

{0\Qix)\v') =V2f^,M^,. (159) 

By defining in the same way the decay constant /t, we introduce into the sum rule the parameters of the 77' and 
P gluonium and the corresponding value of — 6 ~ 7 GeV^ at which has been optimized. In this way, one 
finds that there is no room to include the l contribution, as: 

A « 0. (160) 

One can weaken the constraint by replacing the QCD continuum effect, i.e. all higher-state effects, by the one 
of the L, which should lead to an overestimate of /,,. In this way, one can deduce the upper bound (Fig. 19): 

/, < 16 MeV. (161) 




One can compare the previous result with the one obtained from J/iJ; radiative decays, which gives [pl|: 



Bi,ri' 



r( j/^ ^ 
r( j/v- ^ 7,7') 



{0\Q{x)\l) 



{0\Q{x)\fi') 



(162) 



where the matrix element is controlled by the decay constant of the corresponding particle. Using the experi- 
mental branching ratio, where we take the one for l KKt:, we deduce: 



A ~ 0.23/^. ~ 7 MeV, 



(163) 



in agreement with our findings. Our analysis indicates that the E / i couples more weakly to the gluonic current 
than the 77' (/,,' ~ 30 MeV), and is thus likely to be the radial excitation of the 77' |^. 

The E / i has been interpreted as a bound state of hght gluinos [sot . However, we can consider that the results in do not 

tvour this interpretation, as the inclusion of t he light gluino loops leads to an overestimate of the value of ag{Mz) from r-decays 
3|-[p5|| compared with the world average Wa], This conclusion being confirmed from a global fit analysis of the Z and t hadronic 
decays [p^, and fmm the new ALEPH lower limit of the gluino mass of about 6.3 Gev (95% CL) from the running of Qs and 
four-jet variables |p3]. We also learn from U. Gastaldi that the pp data on KKn can be fitted by one resonance, though the quality 
of the fit is obviously much better within a two-resonance fit. 
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9.3 Meson-gluonium mixing and P lliP'l decays 

Following 1^, we obtain from the evaluation of the off-diagonal two-point correlator, the quarkonium-gluonium 
mixing angle [|[ ^ 

Op ~ 12°, (164) 
from which one can deduce the decay widths of the pseudoscalar giuonium G: 

r(P^77) ~ *^^^^^(^M^j r(r;' ->77) « (1.3±0.1) keV 

r(P^P7) ~ tan^ 6lp^^^ r(?/ P7) « (0.3 ± 0.1) MeV, (165) 

where fej is the momentum of the particle i. We have used r(ry' 77) ~ 4.3 keV and r(ry' p^) ~ (72 ~ 
21) keV. Measurements of the P widths can test the amount of glue inside the P-meson. 

10 Conclusions 

• We have updated the QCD spectral sum rule (QSSR) analysis for computing the masses and decay constants 
of the scalar (section 4), tensor (section 8) and pseudoscalar (section 9) gluonia, using the present values of the 
QCD parameters and the recent progresses in the QCD evaluation of the gluonic correlators. The results of the 
analysis are summarized in Table 1: 



JP(^ Name Mass [GeV] fa [MeV] [GeV] 

Estimate Upper Bound 



0++ 


G 


1.5±0.2 


2.16 ±0.22 


390 ± 145 


2.0 






1.00 (input) 




1000 








1.37 (input) 




600 






3G 


3.1 




62 




2++ 


T 


2.0±0.1 


2.7 ±0.4 


80± 14 


2.2 


0-+ 


P 


2.05 ±0.19 


2.34 ±0.42 


8-17 


2.2 




Eji 


1.44 (input) 




7 : J /ip ^ 





Table 1: Unmixed gluonia masses and couplings from QSSR. 



Our results satisfy the mass hierarchy Ms < Mp « Mt, which suggests that the scalar meson is the lightest 
giuonium state. However, the consistency of the different sum rules in the scalar sector requires the existence 
of a low mass and broad cr-meson coupled strongly both to gluons and to pairs of Goldstone bosons, whose 
effects can be missed in a one-resonance parametrization of the spectral function, and in the present lattice 
calculations. One should also notice that the values of which are approximately the mass of the next radial 
excitations, indicate that the mass-splitting between the ground state and the radial excitations is relatively 
much smaller (30%) than in the case of ordinary hadrons (about 70% for the p meson), such that one can expect 
rich gluonia spectra in the vicinity of 2-2.2 GeV, in addition to the ones of the lowest ground states. 

• We have also computed the masses and decay constants of the scalar quarkonia (section 6). 

• We have used some low-energy theorems (LET) and/or three-point function sum rules in order to predict 
some decay widths of the bare unmixed states (sections 5 and 6), where the results are summarized in Table 2 

g 

^'^We have used the fact that the G couphngs to tttt and KK are negligible as indicated by the GAMS data M 
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Name 


Mass [GeV] 


TT+TT- [GcV] 


K+K- [MeV] 


r]T] [MeV] 


riri' [MeV] 


(47r)s [MeV] 


77 [keV] 


Ob 


0.75 ^ 1.0 


0.2 - 0.5 


5[/(3) 


S'[/(3) 






0.2 - 0.3 




(input) 














'-' B 


1.37 


0.5 ^ 1.3 


SU{3) 


SU{3) 




43 316 


0.7 ~ 1.0 




(input) 










(exp) 




G 


1.5 


w 


« 


1.1 - 2.2 


5-10 


60 - 138 


1.0 ±0.8 


S2 


1. 


0.12 


5J7(3) 








0.67 


S'2 


1.3 w tt' 


0.30 ±0.15 


5L/(3) 


SU{3) 






4±2 


S3 


1.474 ±0.044 




73 ±27 


15 ±6 






0.4 ±0.04 


S'3 


« 1.7 




112 ±50 


SU{3) 






1.1 ±0.5 



Table 2: Unmixed scalar gluonia and quarkonia decays 



• We have discussed some maximal quarkonium-giuonium mixing schemes, in an attempt to explain the complex 
structure and decays of the observed scalar mesons: 

Below 1 GeV: 

We find that, a maximal mixing between two near-by quarkonium S2{uu + dd) and gluonium ag states around 1 
GeV, can explain the large width of the f7(0.75 — 1.), the narrownness of the /o(0.98) and its strong coupling to 
KK (this latter property enables its production from the radiative decay). This scheme being a QCD-based 
alternative to the four-quark and KK molecule scenarios. 
Above 1 GeV: 

— The /o(1.37) is a superposition of two states, the S'2 radial excitation of the S'2 = {uu + dd) quarkonium state 
and the /o(1.37) coming from a maximal mixing between the radial excitation a'^ of broad low-mass a with the 
quarkonium S'3 and gluonium G. 

— The /o(1.5) satisfying the properties observed by the Crystal Barrel collaboration |lj] (namely large widths 
into 47r, 27r and rjrj'), comes also from a maximal mixing between the radial excitation a'g of a broad low-mass a 
with the quarkonium S'3 and gluonium G (orthogonal partner of the /o(1.37)). Our approach also suggests that 
the G(1.6) seen earlier by the GAMS collaboration is likely an "almost" pure gluonium state as emphasized 
earlier in NV. 

— The /,/(1.7) (if its spin is confirmed to be zero) can be identified with the S'3 radial excitation of the S3(ss) 
state. The dip found in the KK mass distribution by the Crystal Barrel and pp annihilation at rest |6^, |67t 
around (1.5 ~ 1.6) GeV can result from a destructive mixing between the /j(1.71) and its orthogonal partner, 
which is very wide. 

• In the tensor sector, using a QSSR evaluation of the off-diagonal quark-gluon two-point correlator, one 
finds, that the quarkonium-gluonium-mass mixing angle is small, of the order of 10° [^, which can exclude the 
identification of the /j(1.71) as a 2++ gluonium Q, but favours the gluonium nature of the observed C(2-2), where 
the total width satisfies our upper bound (section 8). However, due to the small value of the QCD continuum 
threshold, which is about the mass squared of the radial excitation, we expect to have a rich population of 2++ 
gluonia in this 2 GeV region. Our result should stimulate a further search of these states and a test of the 
(non-)existence of the gx seen earlier by the BNL collaboration in the OZI suppressed 7r~p (fxpn reaction, 
but not observed in previous J/^ radiative decay data [76| . 

• In the pseudoscalar sector, the quarkonium-giuonium mass-mixing angle is also small (about 12° [^), which 
combined with the decay widths of the ry'(0.96) allows us to predict the gluonium into 77 width and radiative 
decays (section 9). Finally, we found that the E/ l{1A4) is weakly coupled to the gluonic current, which can 
favour its interpretation as the radial excitation of the ?7'(0.96). 



A further improvement of our mass prediction is under way. 
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